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ABSTRACT 

This  report  describes  a  study  conducted  to  estimate  the  torques  acting 
on  electrically  conducting  shells  rotating  in  a  magnetic  field.  The  basic 
electromagnetic  expressions  leading  to  an  evaluation  of  torque  are  given  using 
various  approaches  with  comparisons  for  the  spherical  shell.  Application  of 
these  techniques  to  cylindrical  shells  with  flat  and  hemisphere  shell  ends  is 
escribed  The  torques  both  slow  and  alter  the  direction  of  the  rotating  mo¬ 
tion,  and  the  resultant  decay  time  constants  can  be  in  the  order  of  days. 
Factors  affecting  this  decay  time  and  influencing  shell  design  are  examined. 
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INTRODUCTION 


HISTORICAL  SUMMARY 

After  the  discovery  of  electromagnetic  induction  by  M.  Faraday  and 

J.  Henry  (1S41)  the  slowing  down  effect,  which  conductors  experience  when 

[l]  * 

rotating  in  a  magnetic  field,  was  studied  by  H.  Hertz  (1896).  Shortly 

r  2  3  ] 

after,  additional  analytic  treatment  was  described  by  R.  Gans  (1903). L  ’  1 

These  results  for  spherical  shells  also  appear  in  the  relatively  modern 
f  4l 

text  by  Smythe. L  1 

Effort  during  the  first  half  of  this  century  concentrated  on  generator 
and  motor  design.  While  studies  in  electromagnetism,  terrestial  magnetism, 
and  mechanics  appeared  in  various  texts,  a  rekindling  of  interest  in  the 
basic  questions  concerning  the  interaction  of  conducting  bodies  moving  in 
magnetic  fields  awaited  the  arrival  of  the  space  age  around  1955.  Then  a 

r5  6  7] 

number  of  investigations,  for  example,  Rosenstock,  Vinti,  and  ZunovL  ’  ’ 
considered  the  problem  in  terms  of  a  spinning  satellite  moving  in  an  earth 
orbit.  Yet  until  the  present,  very  little  was  done  analytically  on  other  than 

a  few  regularly  shaped  bodies  such  as  spheres,  cylinders,  and  prolate 

,  ..  [8,  9,  10] 

spheroids. 

A  good  deal  of  effort  has  been  given  to  using  observed  satellite  data  to 

arrive  empirically  at  decay  and  motion  perturbation  estimates.^1,  12  ^  In 

[13] 

this  regard,  a  rather  comprehensive  treatment  recently  appeared  by  Yu. 


Numbers  in  brackets  refer  to  References  cited  at  the  end  of  the  report. 


L 


Finally,  a  relatively  large  effort  also  has  been  given  to  utilizing  the 

interaction  of  the  magnetosphere  with  on-board  magnetic  sources  in  order 

[ 14  15] 

to  effect  attitude  control/  ’  In  dealing  with  objects  in  orbit,  the  inter¬ 
action  of  the  magnetosphere  need  not  only  be  with  the  induced  eddy  currents. 
It  may  also  react  with  existing  current  or  magnetic  sources  on  the  object  as 
well  as  with  the  magnetic  field  induced.  There  may  also  be  reactions  of 
ambient  electric  fields  with  induced  electric  fields. 

THE  SCOPE  OF  THE  PROBLEM 

When  a  conductor  and  a  magnetic  field  undergo  relative  motion*,  an 
electromotive  force  is  developed  in  the  conductor.  If  closed  conducting 
paths  are  present  through  which  a  change  of  magnetic  flux  occurs,  then 
eddy  currents  flow  in  the  conductor.  The  further  interaction  of  this  current 
flow  with  the  magnetic  field  produces  torques  which  in  general  both  dissipate 
the  rotational  energy  and  cause  an  orientation  change  in  the  rotation  axis. 

In  addition,  the  currents  if  sufficiently  strong  produce  a  magnetic  field  which 
can  appreciably  perturb  the  background  magnetic  field. 

The  analytic  effort  of  interest  here  is  that  concerned  with  deriving 
expressions  for  the  resulting  torques,  and  also  with  the  decay  time  con¬ 
stants  which  result,  rather  than  the  over-all  motion  problem  to  account  for 
all  magnetic  torque  components  and  the  resulting  motion  about  the  center  of 
mass  as  the  object  moves  in  orbit. 

It  will  also  be  our  concern  here  to  deal  with  modified  cylinder  shapes 
and  to  discuss  the  design  problem  for  realizing  desired  motion  lifetimes. 


♦This  includes  the  case  of  a  stationary  conductor  and  a  magnetic  field  varying 
in  time-space. 
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SECTION  n 


ANALYTICAL  APPROACHES 

BASIC  ANALYTICAL  MODEL 

To  provide  a  background  for  the  methods  of  calculating  eddy  current 
torques,  we  begin  this  section  by  describing  the  basic  electromagnetic  and 
mechanics  equations  appropriate  to  establishing  mathematical  models.  The 
following  symbols  apply  to  the  equations  given  in  the  report: 


Electromagnetic  s 
6  =  electromotive  force 
E  =  electric  intensity  vector 
H  =  magnetic  intensity  vector 
D  =  electric  flux  vector 
B  =  magnetic  flux  vector 
$  =  flux  linkage 
p  =  charge  density 
F  =  current  density  vector 
a  =  conductivity 
€  =  emissivity 
p  =  magnetic  permeability 
c  =  speed  of  light 


Mechanics  -  Geometry 
M  =  mass 

I  =  moment  of  inertia 

to  =  angular  velocity  vector 

V  =  velocity  vector 

r  =  distance  from  center  of  mass 

f  =  force  vector 

T  =  torque  vector 

£  =  path  length  on  object 


Using  a  Gaussian  system  of  units,  the  applicable  basic  equations  are: 


V  x  E  =  -  - 


1  0_B_ 
c  3t 


a  E  +  -£  (V  x  H)  J  +  p  V 


3 


(3) 


V  x  H 


r  + 1 iv. 

c  c  d  i 


V  =  w  x  r 


(4) 


Equations  (1)  to  (3)  apply  whether  measured  in  an  inertial  frame  at  rest  with 
the  laboratory,  or  fixed  with  respect  to  an  inertial  observer  instantaneously 
at  rest  with  respect  to  a  point  in  the  moving  object.  For  an  inertial  frame 
at  rest  for  some  point  on  the  object,  the  following  constitutive  equations 
apply:  ' 


D'  =  e  E' 

(5) 

B'  =  pH' 

(6) 

Using  unprimed  quantities  to  apply  for  a  fixed  frame  and  taking 

2 

H  ^  1,  €  =  1  (for  a  metal)  and  (V/c)  «1  we  get,  inside  the  metal,  for 

[6] 

our  purposes 


Also  for  (V/c)  «  1 


D  s  E 
B  =  H 

i  =  i1  +  p  V 
E*  =  E  +  (V  x  H)  p/c 


(7) 

(8) 

(9) 

(10) 


In  addition,  the  contribution  of  the  convection  currents  p  V  in  regard 

to  producing  torque  can  be  neglected.  Their  elimination  in  this  respect  can 

be  justified  both  due  to  the  relatively  small  contribution  to  i  and  also  a 

possible  cancellation  for  torque  due  to  geometrical  symmetry.  The  only 

[  6] 

effect,  which  is  very  small,  is  a  turning  torque.  J  However,  in  general 
we  take 


4 


1 


(11) 


=  f«  +  pv 

V  =  a  [l  +  £  (V  x  H)]  (12) 

In  respect  to  the  discussion  in  the  introduction,  it  is  interesting  to 
recall  the  two  mechanisms  for  inducing  an  electromotive  force  6  about  a 
conducting  path  S.  Using  n  =  1  we  have 


=  '  5  [(h)  v  =  0  +  H  B  =  cons,imt  ] 

1  fit. 

c  dt  (13) 

where  d$/dt  is  the  rate  of  change  of  flux  linking  path  S. 

We  are  concerned  here  with  (3$/3t)  B  =  constant  due  to  motion  as 
well  as  (3$/3t)  V  =  0  if  3H/3t  ^  0,  that  is,  if  the  background  magnetic 
field  is  disturbed. 

If  the  background  magnetic  field  is  fixed  and  the  eddy  currents  do  not 
appreciatively  alter  it,  then  we  may  take  3H/3t  =  0  in  Equation  (1).  This 
allows  V  x  E  =  0  which  results  in  a  convenient  approach  for  the  spherical 
shell  described  under  "The  Spherical  Shell.  " 

The  angular  motion  w  can  be  resolved  as 

w  =  wti  +  wi  (14) 
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where 


wjj  is  parallel  to  H 
cdj  is  normal  to  H 

The  effect  of  coj  j  is  to  cause  charge  separation  toward  the  conductor 
boundary,  that  is  radially  from  the  rotation  axis,  resulting  in  an  electric 
field  E  which  balances  out  the  originally  induced  electric  field  caused  by 
the  motion,  hi  other  words  with  3H/3t  =  0  and  taking  Vi  i  =  un  x  r  then 


E'  =  E  +  -  (  V„  x  H  )  =  0 


(15) 


Using  Equation  (12)  then  gives 


i*  =  0 


(16) 


With  respect  to  the  fixed  frame  the  zero  torque  result  also  obtains  as 
can  be  seen  from  the  force  expression  using  Equation  (11)  for  i  as  follows 

f  =  pE  +  ixpH  =  pE  +  [pV  +  cr(E  +  Vxp  H)]  x  p  H 

=  p(E  +  VxpH)  +  a(E  +  VxpH)  =  0  (17) 

from  Equation  (15). 

This  is  equivalent  to  saying  that  under  Wjj  action,  no  closed  i  loops 

occur.  Here  V  ‘  i  /  0  and  amounts  to  3  D/3  t  /  0.  The  effect  is 
referred  to  as  electrostatic  shielding.  For  j  the  charge  separation 

dissipates  rapidly  along  closed  current  loops.  This  in  essence  amounts  to 
assuming  3  D/3  t  =  0  ,  since  using 


D  =  4  7T  p 


(18) 


6 


implies  no  charge  build-up.  More  exactly,  from  Equations  (3)  and  (18)  with 
D*  ,  we  have 

„  .  1  av  •  D'  8  p 

V  ’  1  "  "  4tt  8 1  8 1  (19) 

The  rates  8  p'/ 8  t  are  electromagnetic  and  so  much  more  rapid  than 
the  motion  rate  u  .  Thus,  taking  for  the  transient  i'  =  a  E',  Equation  (18) 
allows  simply 


8  p'  iirap' 

rr  e  (20) 


that  is 


P' 


e/4  ra 


(21) 


—18 

For  example,  t  for  aluminum  s  0.2  x  10  sec.  Thus,  for  a  conductor 
such  as  aluminum,  the  time  constant  t  is  extremely  short. 


Hence,  any  tendency  for  charge  build-up  is  rapidly  dissipated  and  in 
the;  steady  state  (w  ~  constant)  causes  no  opposing  electric  field  as  in 
Equation  (16).  Notice  from  Equation  (1)  that  when  8H/8  t  /  0  ,  an  electric 
field  can  be  developed  even  though  no  charge  build-up  occurs.  Thus  we  may 
summarize  by  stating  that  an  electric  field  is  produced  to  oppose  that  caused 
by  the  motion  (V  x  H/c)  if 

(1)  i  causes  8  H/8  t  /  0  (i.  e. ,  changes  in  the  total  magnetic  field). 

(2)  i  causes  a  build-up  of  charge. 


7 


m 


x 


In  the  case  of  only  wjj  present,  case  (2), as  we  have  noted,  results  in 
i*  —  0  and  so  no  torque  is  produced.  *  In  such  a  case  the  opposing  electric 
field  is  generated  by  the  charge  build-up. 


Considering  Equations  (3)  and  (11)  leads  to 


V  x  H  =  —  (e  + 
c  ' 

We  see  that  3  D/3 1  =  0  requires 


-  •  v  x  H  ) 
a  / 


3D 
c3  t 


,  Tr 

+  —  p  V 
c 


(22) 


1*2.  E 
c 


» 


e  3E 
c  3  t 


or  for  harmonic  time  variation 


(23) 


oVf  »  £ 

a  condition  satisfied  for  good  conductors  and  frequencies  up  to  1010  me. 

In  this  discussion  the  condition  is  valid  and  results  in  the  so-called  quasj- 
stationary  case  with  no  displacement  current. 

Taking  the  case  of  Equation  (22)  we  have 


v2h 


4jro-p  3H  47ro-p  /_  .r.  rr,\ 

-r  ft  -  —  (vx  <v  x  HV 

c  c 


We  consider  E 


Ei + 


such  that  we  associate 


(24) 


(25) 


♦Unless  wjj  is  strictly  constant,  current  transients  occur.  The  torques 

developed  from  the  resulting  field  disturbances  are  negligible,  however, 
with  respect  to  the  torque  due  to  the  steady  state  co  j  motion. 


8 


and 


~r~  V2  H  =  -  apV  x  E  (26) 

4  7T  2 

For  example,  an  aluminum  object  of  characteristic  length  i  =*  125  cm  and 
with  V  on  its  surface  between  0  and  about  12  cm/sec  (corresponding  to 
<y  £*  2  jt/30  rad/sec),  has  V  x  (V  x  H/c)  ~  VH/c£  and 
c/Axap.  V2H  ~  cH/4ira/i£  with  equal  order  effects  on  dll/dt.  Then  E^ 

and  E„  should  be  comparable  where  we  have  associated  E  with  the  motion 

—  —  —  2—  ^ 

V  x  H/c  and  E  with  V  H,  that  is  the  effect  of  the  presence  of  the  body 

Ld 

in  disturbing  the  background  magnetic  field. 

We  shall  conclude  this  discussion  by  developing  the  basic  torque 
expressions. * 

From  Equations  (4)  and  (11), 


and 


(27) 


(28) 


where  the  differential  torque  dT  is  developed  from  the  force  occurring 
on  a  differential  volume.  Then  with  (i  =  1 


*By  neglecting  the  p  ?  term  in  I  we  are  in  essence  taking  E  to  be  effec¬ 
tively  E  in  E  =  E  +  E„  where  E  is  due  to  the  component  of  co 
J  2  _  1  2  1 

parallel  to  H,  i.  e. ,  the  component  due  to  charge  build-up. 
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s 

T  =  J  dT  =  l  frxjo-  -J  x  H  +  eJ  x  h|  dV 

Vol  Vol  \  '  (29) 

As  noted  previously,  for  w  =  wjj  .  T  =  0  (no  flux  linkage).  Then  expanding 
(w  x  r)  x  H  and  using  Wjj  x  H  =  0,  we  get 


l  I 


r  x  <  (fa>  ‘cH)  r  +  E  x  H  dV  =  0 


Since  r  can  have  arbitrary  orientation  with  respect  to  H  .* 


=  w  •  H  \  -  ,  - 

E  =  - 1 — -  I  r  ;  for  w  =  ui 


which  is  an  alternate  form  of  Equation  (15). 


H  =  H|  +  Hjj 


where 


Hjj  is  the  component  of  H  parallel  to  w 
H  |  is  the  component  of  H  normal  to  w 


Similarly  let 


E  —  Ei  +  E 


1  +  "II 


It  is  to  be  recalled  (see  Equation  (17))  thatjhe  E  of  Equation  (31)  is  due  to 
charge  build-up.  We  have  neglected  the  p  V  term  in  the  torque  expression 
here  since  it  does  not  effectively  contribute. 
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where 


i 


v  _  w  •  H  :  _  wjHnlr 

I!  ‘ -\s—  r  - '  - 


from  Equation  (31)  using  1 1  to  stand  for  magnitude  and  let 


Ei  =  E  '  Eil 


(34) 


(35) 


We  note  if  co  =  ,  then  (o>  •  H)  =  0  so  that  E  =  E 


Then, 


’si 


T  =  -  J  r  X 

Vol 


E,  - 


1 


(H  *  r)  (i) 


1 


x  H  |  dV 


3  f 


r  x 


Vol 


We  can  write  Equation  (36)  as 


r_  -i 

<7  P  - 

i — 

X 

a 

i _ 

dV  g  Jrx( 

C  T 

(H  •  r)  co  x  H 


(36) 


dV 


where 


5  ’l  y-r  x  (Ej^x  H)dV  ♦  2  f0 
Vol 


T  =  (&)  x  H)  x  (H  •  r)rdV 

°  2c  vi 

Vol 


(37) 


(38) 


is  an  expression  for  the  complete  torque  on  a  sphere  shape  if  V  x  E  =  0 
is  assumed  as  well  as  9D/8 1  =  0. 


11 


The  Spherical  Shell 


We  begin  the  discussion  using  an  approximate  but  convenient  method 

[  6l 

based  on  Vinti.L  This  method  also  gives  insight  into  torques  for  other 
shapes  discussed  in  Section  in. 

The  assumption  made  is  that  V  x  E  =  0.  This  amounts  to  saying 
that  the  magnetic  field  is  not  appreciably  perturbed  by  the  rotating  conducting 
sphere  shell.  We  define  the  inner  radius  of  the  shell  by  a  and  the  outer 
radius  by  b,  and  call  the  thickness  h  =  b  -  a.  As  noted  previously  the 
torque  is  given  by 


T=  J[rx 
Vol 

♦with  V  x  E  =  0  (i  e  8H/9t  =  0)  then 

V  x  i  =  -  V  x  (V  x  H)  . 

For  H  constant  and  uniform  and  o>  constant,  we  find  expanding  the  curl 
of  the  cross  product. 


(I  x  H) 


dV 

c 


(39) 


Then 


Vx(VxH)  =  -Hxcu 


V  x  i  =  ^  (a)  x  H) 


(41) 


(42) 


or 

i  =  ~  (w  x  H)  x  r  {■  V  <p  (43) 


♦The  implication  is  that  E  =  W  with  a  scalar  function,  since 
V  x  (V*)  =  0. 
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We  note  that  V  •  i  =  0  using  3D/3t  =  0,  and  with  Equation  (43) 
and  V  •  [  (w  x  H)  x  ?]  =  0  ,  then 

2 

V  <p  =  0  (44) 

that  is,  4>  satisfies  Laplace's  equation. 

Considering  the  surface  boundary  conditions  due  to  continuity  of  the 
normal  i  component 

3  (b  i 

—  =  S7  <p  \  gurface  =  0  and  the  geometrical  radial  symmetry 
applicable  with  the  sphere,  there  results  V<p  =  0  or 


1  =  2c  (w  x  H)  x  r 


(45) 


We  then  obtain  from  Equation  (39)  by  expanding  the  triple  product  and 
simplifying  results  in  a  single  term, 


T  =  Tq  =  —  (w  x  H)  x 
2c 


I  H  *  y  r  r  dV  j 
L  Vol  J 


(46) 


Now  for  shell  sphere  taking,  *  r  =  x  i  +  x  i  +  x  i 

X  JL  Lt  Ct  o  o 


* 

The  i  ;  i  -  1,  2,  3  are  unit  vectors  in  the  3  coordinate  directions. 

Equations  (47)  come  from  \  x.x.  dV  —  0,  i  /  j  and 

J  1  J 

Vol 


C  x2  dV  =  \  f  r2  dV 
J  i  3  J 

Vol  a 


1,  2,  3. 
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/ 


I 

Vol 


rrdV  = 


=  %  (■ b5  - 
15  \ 


^  A  A  A  A  A 

11  +  ^2  *2  *3  *3 


i 

A 


With  H  -  H,.,  ♦  H2t2  +  H3f3, 


H  '  I  ffdV  =  IS  (b5  -  v)  H 


Vol 


Thus 


T  = 
o 


2tt  cr  /  5  5  ,  TT. 

15  -2  \ b  -  a  I  <"  x  «) 

c  ' 


5)(ii 


x  H  =  M  x  H 


where  we  may  take 

M  =  magnetic  moment  of  the  sphere  shell. 
We  define  the  constant  factor 


,  27 ro  /  ,5  5  \  la 

k  =  - *  I  b  -  a  =  - ; 


15c 


4pc 

27rffa'Si 


3c 


(47) 


(48) 


(49) 


where  I  is  the  moment  of  inertia  of  sphere  shell,  p  =  uniform  mass  density 
and*  for  h/a  «  1  , 


In  general  for  h/a  «  1,  bn  -  a11  =  a“  A  (nh)  or  b“  A  (nh). 


n  -  1 


n  -  1 
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=  5b4h 


(50) 


Thus 

To  =  k[(S  x  H>  x  h]  (51 

We  may  resolve  Tq  into  components  parallel  and  orthogonal  to  w  denoted 
respectively  as  T  ^  and  (D  for  drag,  T  for  turning). 

Then 


Substituting  from  Equation  (51)  and  expanding  the  triple  cross  product 


foD  ’  -“H2  -(fu  '  B/_ 


(v  =  -  k  j  Hj^  J  ^  w 


(53) 


and 


oT 


-  k“  (i  ■  H)[i  x  (B  x  fo.) J  k“  lH|ll  H1 


(541 


with  Equation  (48)  it  is  to  be  noted  that  Equation  (51)  becomes 

2 


oD 


2tt  |h  j  |  (s  (b5  *  a5) 

- - - oj 


15c 


(55) 


Equation  (55)  is  new  described  using  an  alternate  approach  which  is  less 
restrictive  and  applicable  to  other  object  shapes. 

Let  us  consider  the  following  geometry  for  the  solid  sphere  of  radius  b. 
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Let  {x,  y,  z  }  be  a  fixed- in-space  coordinate  system. 

{|,  h,  Z}  be  a  fixed-in-the-body  coordinate  system. 


We  now  reformulate  the  model  of  a  spinning  shell  in  a  fixed  field  to  a  shell 
in  a  field  which  changes  in  time  to  correspond  to  the  relative  motion,  since 
torque  depends  only  on  relative  motion.  We  note  also  that  even  though  H 
is  taken  normal  to  co  in  the  diagram  so  as  to  get  only  a  damping  torque,  the 
model  allows  a  general  orientation  between  H  and  o>  as  in  the  previous 
calculation  for  torque.  For  the  motion  assumed  we  can  take  Z  =  z.  Then 

H„  =  (H  cos  cutlf. 

£  \  x  /  £ 

H  =  -  ( H  sin  w  t)  1 
V  \  x  /  n 


Hz  =  0  (56) 

so  that  H  =  H  +  H 
€  1 
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In  complex  form  these  become 


H  =  H  e  w  ’  H  =  Re  H 

£  x  £  I 

-tot . 

H  =  -  iH  e  H  =  Re  H 

r\  x  t]  i) 


The  method  now  uses  the  formulation  of  the  magnetic  moment  based 

f  16l 

on  Landau  and  Tifshitz.  L  J 

Let  M  =  total  magnetic  moment  acquired  by  a  conductor  in  a  mag¬ 
netic  field.  In  general,  the  torque  on  the  conductor  in  the  magnetic  field 
H  is  then 


T  =  M  x  H  (58) 

which  compares  to  Equation  (48)  for  the  sphere.  We  take 

3  3 

M  =  /  M.  i.  ,  H  =  V  H.  l. 

l-j  xi  U  J  J 

i  =  1  j  =  1  (59) 


Then  M.  can  be  written  as 
i 

M.  =  Vo:.,  (w)  H, 
i  lk  k 

where 


(60) 


V  =  volume  of  the  conductor 

w  =  angular  rate  of  conductor  with  respect  to  field  direction 


Va.^cu)  =  magnetic  polarizability,  a  tensor,  where  a 
is  symmetric  and  depends  on  the  body  shape 
and  orientation  of  the  external  field  but  not  on 
the  volume. 

For  the  sphere  (solid)  a  is  a  scalar*  a  6^  and  we  have 


M  =  V  a  H 


(61) 


with  a  =  a  +  ia  For  the  sphere,  the  same  a  applies  whether  the 

X  c* 

H  field  is  Hj^  or  Hu  so  that 


I! 


M  ,  =  VH  (o:,  cos  cut  + 

£  X 


^  cos  cut  +  a2  sin  cut^ 


M  =  V  H  / -a,  sin  cut  +  a  cos  cut } 

t\  x  ^  1  2  / 


(62) 


M  =  0 
z 


Using 


M  =  M  cos  cut  -  M  sin  cut 
x  £  V 


M  =  M,  sin  cut  +  M  coscut  (63) 

y  I  V 


*6.,  is  the  Kronecker  deita.  For  the  sphere  the  polarizability  per  unit 
lk 


volume  is  a  =  a  t-  ia  = 


where 


k  = 


1  +  i 


-  — - -  and  6  = 


.  -  JL  [i  . 

8  7T 

is  the  skin  depth.  It  is  to  be  noted  also 


—  +  —  cotbk 
,2. 2  bk 
b  k 


(2iracu)  l/2 


that  since  a  is  complex,  a  retardation  exists  between  M  and  H  and  so 
non-colinearity  and  non-zero  torque. 
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We  get 


M  =  VH  a 
x  xl 


M 

y 


V  H  a 
x  2 


M  =  0 
z 


(64) 


Then  with  H  =  H  i  and  Equation  (61), 

X  X 

T  =  0 
x 


T 

y 


T 

z 


0 


-M  H 

y  x 


-  V  H2  a 
x  2 


(65) 


If  the  skin  depth*  o  »  h  then  a  =  ^  so  that  with  V  =  ^  Trb2, 

Lu  U  O 

10c 

.  ,3  2  2  *  TT  2  5- 

4irb  H  b  oui  27rH  ob  0) 

T  _ _ x  z  _  x 

Z  30c2  15c2  (66) 

To  approximate  the  thin  spherical  shell  when  h/a  «  1,  we  subtract 

the  result  for  T  at  a  ,  the  inner  radius,  from  the  result  T  at  b,  the 
z  z 

outer  radius.  This  gives  for  the  spherical  shell  from  Equation  (66), 


T  =  - 
z 


H 


/,  5  5\  - 

alb  -ala) 


(67) 


For  a  shell  5  »  h  (the  shell  thickness). 
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Equation  (67)  agrees  exactly  with  Equation  (49)  since  in  our  case 

H  =  H,  =  H  r  . 

xx 

The  form  given  by  Equation  (67)  thus  agrees  with  Equation  (55)  which 
was  derived  neglecting  the  effect  of  9  JI/9  t  /  0  ,  if  present. 

For  a  spherical  shell  thick  enough  and  with  the  H  geometry,  the 

assumption  9  H/9  t  =  0  is  reasonable.  The  formula  for  a  very  thin  shell 
[4] 

given  in  Smyth  includes  the  9  H/9  t  effect  and  restricts  H  normal  to  co. 
Converting  this  result  to  our  terminology  gives. 


T  = 


C  co 

1 


1  +  C2w 


where 


Cx  =  (2/3)7?M2H2ab4  h 

2  2  2,2 

_  M  b  cr  h 

2  9 


(68) 


For  us  (i  =  1  and  taking  9  H/9  t  -  0  is  to  take  C  =  0.  Using 

u4  _  b5  -  a5  .  2ttH2(t  (b5  -  a5)  ,  .  , 

b  =  — ^ -  then  gives  C,  =  - jg -  .  Applying  this  to 

the  Gaussian  system  of  units  then  produces  agreement  with  Equations  (55) 

and  (67). 


We  close  this  section  by  noting  that  a  recent  investigation  by  Halverson 
r  g  j 

and  Cohen1  Jhas  included  the  generalization  of  considering  both  rotation 
with  co  and  B  arbitrarily  related  as  in  Equations  (52)  or  (59)  as  well  as 
9H/  /  0  as  in  Equation  (68). 

The  model  of  using  the  assumed  equivalence  of  a  time-varying  f'.eld  and 
fixed  conductor  for  the  moving  conductor  in  a  fixed  field  is  applied. 
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Using  the  vector  potential  A  (see  References  (4)  and  (8))  where 

V  x  A  =  B  (od) 

then 

r72  7  -  4^a  dA 

c  d  t  (70) 

and 

t-  co  dA 

1  “  "  4t  dT  (71) 

We  note  Equation  (70)  also  applies  to  E,  H  and  i  . 

The  torque  given  by  Equation  (29)  or  (39)  can  be  written  in  the  equivalent 

form 


T  =  -  - 
c 


j  ?x[^x  (v  x  a) 
Vol 


dV 


(72) 


Solving  Equation  (70)  with  the  appropriate  boundary  conditions  for  the 

sphere  shell  results  in  a  general  formula  with  rather  complicated  coeffi- 

[8] 

cients.  For  a  good  conductor  or  with  u>  low  (with  q  —  0,  see 
Reference  (8)]  or  a  good  dielectric  or  high  cu  (q  —  ®°)  ,  the  form  of  the 
coefficients  simplify,  with  w  =  wk  and  £  (cJ,  T3)  =  a  ,  to  give  for 


(73) 
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for  q  >  100  (q  —  «°) 


where 

q  =  b  (2  P)1//2;  P  =  ,  t  =  resistance/unit  length 

T 

p  =■  b/a 

For  a  thin  shell  p  —  1,  so  that  taking  h  =  b  -  a  with  surface 
resistivity  £  =  r/h  remaining  finite,  then  with  q  arbitrary, 


TC"b  2~2  2  [(“  **  +  **  w  bl)  sin  2a  +  £  ic(2sin2a)l 

+  P  w  b  L 


9£  +  P 


For  Bj^  w  (i.  e. ,  a  =  ir/2) 


2-4 
6t  B  £wb 

n  2  ^  2  2.2 
9£  +  p  «  b 


2-4 
6t  B  ocub  h 

n  2  .2  2  2 

9  +  p  b  a  h 


Equation  (67b)  is  the  form  given  by  Equation  (68)  with  £  =  l/oh  . 

The  i  ,  j  ,  and  k  components  of  T  in  Equation  (73)  ,  (74)  or  (75)  are 

2 

-C  sin  2a  ,  C  sin  2a  ,  and  2  C  sin  a  respectively.  These  forms  also 

X  m  X 

apply,  aside  from  the  specific  values  for  and  C 2  ,  for  the  general  case 
with  no  restrictions  on  q  or  h. 
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We  see,  as  before,  that 


(1)  T  =  0  if  or  =  0  (i.  e. ,  o>  and  B  are  parallel) 

(2)  The  i  and  f  components  of  T  cause  o>  to  precess  about  B  . 

The  situation  in  (2)  is  consistent  with  the  simplified  form  of  TqT  in 
the  H|  direction  given  in  Equation  (54).  Except  for  perfect  conductors, 

JL 

the  vector  to  then  precesses  in  the  steady  state  until  it  becomes,  for  fixed 
B  ,  colinear  with  B.  As  ai  approaches  B  ,  energy  is  transferred  to 
ct’jj  which  is  dropping  due  to  damping  of  o>  such  that  coy  remains  constant 
for  constant  B. 


/ 


SECTION  III 


application  to  cylinder-type  shells 


NON-SPHERICAL  SHAPES 


As  noted  the  torque,  in  general,  can  be  given  as 


T  =  - 
c 


1  r  -  s  r 

-  \  r  x<  cr 
c  J  ) 


kjJULl L+e  X  H  /  dV 


where  the  effect  of  convection  currents  can  usually  be  neglected.  Expanding 
the  triple  product  then  gives 

T  - 1  Jr  x  |7l  +  J  X  h]  dV  -  I  Jr  x^i-a) 

Vrvl  *-  \  /  J  L  / 


x  H  dV 


Using  Equations  (33)  and  (34)  then  aUows  for  T,  in  general,  as 


dV  +  2  T 


where 


L  '  ("  *  H)  x  (5  '  I  ”r  4 

\  Voi  / 


In  the  case  of  a  sphere  shell  T  =  Tq  which  implies 


?  j  -r  x  (ej_  x  hW  -  -  T0  . 

tr-1  '  f 


24 


« 


To  say  this  in  general  would  be  to  imply  that  V0  =  0  in 

I=^^wxH)xr  +  V(^  (81) 

which  is  not  true  for  non-spherical  body  geometries. 

Thus,  approaches  using  equations  such  as  Equation  (70)  which  take 
account  of  the  field  and  currents  internal  to  non-spherical  geometries  must 
be  used.  For  cylinder  shapes  T  values  can  result  which  do  not  differ 
greatly  from  Tq  .  In  addition,  the  dyadic  form  ^  r  r  dV  in  Tq  is  useful 
for  comparing  different  cylinder  end  configurations.  (See  Equation  (83)  and 
section  entitled  "Hemisphere  End  Cylinder". ) 

We  also  note  that  if  conditions  (1)  and  (2)  (see  page  9)  do  not  hold,  then 
as  also  noted  in  Reference  (17),  we  may  drop  E  and  take 


-o-/- 
i  =  -  (w  x  r 
c  \ 

To  do  this  would  give  for  H  , 

T  =  2  T  (83) 

o 

which  in  the  sphere  case,  at  least,  is  in  error  by  a  factor  of  2. 

Thus,  although  V  x  E  =  0  is  assumed  from  3H/8t  =  0  to  get 
T  =  Tq  for  the  sphere,  charge  movement  etc.  does  not  allow  dropping  of 
the  E  itself. 


x  H 


(82) 


Finally  when  oj  —  wy  then  T^  —  0  .  With  T  —  0  in  general  for 
such  orientation,  then  we  might  use  the  approximation 


since 


T  =  K  T 

o 

and  E^  are  not  in  general  zero. 


(84) 
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RIGHT  CIRCULAR,  CYLINDERS 


Initially  in  order  to  ’ify  the  field  equations,  a  spin  axis  co  normal 
to  H  is  assumed.  This  a.  u  ms  comparison  with  more  complicated 
analysis.  Here  the  assumptions  of  both  zero  displacement  current  3D/3t  =  0 
and  nonperturbed  field  9H/8t  =  V  x  E  =  0  are  made. 


Considering  the  field  induced  in  a  current  path  moving  with  the  body 
we  write  from  Equation  (10) 


9H 

VXE'  =  *  V  X 


£  (v  x  h) 


(85) 


With  V  x  E  =  n  9  H/9 1  =  0,  V  |  the  component  of  V  normal  to  H,  f 

A 

the  cylinder  length  at  direction  £  ,  and  neglecting  end  effects,  there  results 


E'  =  —  V.  £ 
c  1 


M  H 


r  a)  sin  a  £ 


(86) 


where  a  =  }  (V,  H)  . 

For  a  path  length  f  ,  the  voltage  5  becomes* 


H  „ 

&  =  —  £  r  oj  sin  a 
c 


(87) 


An  approach  similar  to  the  one  outlined  here  can  be  based  on  using 
V  =  Voltage  =  d^/dt ;  <p  =  flux  linkage. 

P  =  power  dissipated  =  -V2/R  ;  R  =  resistance  of  induced  current  paths. 

T  =  torque  =  l/w  dV/dt  =  1/w  P. 

Such  an  approach  can  be  used  to  develop  torque  values  for  tumbling  cylinderJv^ 
However,  such  methods,  although  simple  to  construct,  can  be  in  error  (e.  g. , 
page  30)  since  they  do  not  account  completely  for  the  induced  current  paths  and 
boundary  conditions  as  do  methods  based  on  Maxwell  field  equations. 
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By  calculating  the  path  resistance  next,  a  value  of  ~i  and  then  T  is  obtained 
using 


'if 


rx(ixH)dV 


With  a  and  b  the  inside  and  outside  radii  respectively  for  a  cylinder  of 

[9] 

length  £  ,  Hooper  using  this  approach  gives,  with  conducting  end  plates 
(in  MKS  units), 

,  2  2  2  ( ,  4  4^ 

„  nfi  R  £  wo\d  -ay 

4  (£  +  a  +  b)  —  (88) 

or 

2 „2 . 2  ,3  . 

T=  H  £  uirb  h 

(£  +  a  +  b)  (89) 

using  h  =  b  -  a  and  the  footnote  on  page  14,  and  without  end  plates, 


2  2  2  (,4  4  ^ 

T  =  -  T  H  e  w  ~  a  / 
4[£  +  it / 4  (a  +  b)j 


or  as  in  (89) 


2„2  2 

• r  H  H  £  iiiirb  h 
£  +  t/4  (a  +  b) 


Similar  calculations  on  the  sphere  shell  give 

-  32  2  2  (,  5  5\ 

T  =  -—  p  H  w  cr  \b  -  a  j 


(90) 


(91) 


(92) 


which  agrees  very  closely  with  Equation  (67). 

However,  Smythe's  analysis  (Reference  4)  accounting  for  3H/3t  ^  0 
gives  for  w,  the  spinning  cylinder  shell  without  end  plates,  Hj^  c3  ,  and 
h  small 


27 


T  = 


(93) 


4  ir  w  a  h  b3  B2  i 

.  ,  22,2,22 

4  +  co  n  b  h  or 

Neglecting  9H/9t  reduces  Equation  (93)  to  (in  MKS  units) 

T  =  -  t  u  a  h  b3  B2  t  (94) 

which  agrees  with  (91)  if  the  term  t  /a  (a  +  b)  is  neglected  in  the  denomi¬ 
nator.  Since  (94)  agrees  with  a  third  calculation  given  below,  it  appears  that 
the  calculation  of  the  resistance  path  in  Reference  9  may  be  in  error. 

We  now  proceed  to  a  calculation  similar  to  that  used  on  the  sphere 
starting  with  Equation  (56).  This  will  be  done  for  a  tumbling  (as  well  as 
spinning)  cylinder.  In  both  cases  o,^K  still  applies. 

We  view  the  geometry  as  shown 


y 


For  a  solid  cylinder  the  magnetic  moment  per  unit  length  is 

M  =  V  a  K 

where  the  tensor  a.,  reduces  to'the  scalar  ad., 
lk  lk 

a  =  a,  +  ia 

x  2 
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Equations  (56)  through  (60)  still  apply.  Now,  however,  we  deal  with  two 
values  of  a,  say  a *  and  a",  where  a'  applies  when  H  is  orthogonal 
to  the  cylinder  axis  and  a  "  when  H  is  parallel  to  this  axis. 

Following  as  before,  we  now  obtain 


=  VRe 


(0^)=  VRe  (o^  +  i  «2')(Hxe 


-icut 


=  VH  (  a:  *  cos  cut  +  of  *  sin  wt 
x  \  1  2  / 


(95) 


M  =  VRe 
V 


=  VRe 


=  VH 


VH 


(V  +  iV)  (-iI/'t) 

^-a^"  sincut  +  ar2"cos  cut) 

(-*1*  sin  wt  +  <Yc°s  «t)  (96) 


Mz  =  0^Hz  =  0^ 


(97) 


using,  in  the  equation,  the  fact  that  a"  =  a '/2  (see  Reference  16). 


With  Equation  (63)  we  have 


M  =  VH 
x  x 


M  =  VH 
y  x 


V  1 


.  2  , 
sm  cut 


\  ,  si 

J  +  a  2  ~ 


sin  wt  cos  cut 


a. 


sin  cu  t  cos  cu  t 


,  /  2  cut 

+  “2  1  "  C°S  X, 


(98) 


(99) 


M  =  0 
z 


(100) 
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T  = 
z 


M  H 

y  X 


,4  : 
rb  a 

7 

6c" 


2  2 

co  sin  co  t  cos  <o  t 


2  /  2 
b  aw  (  ,  cos  co  t 


VH  2  b2  a  co  \r  .2  /  2 

x  irb  era;  A  /,  cos  wt 

- o -  / - o~  sm  a>t  cos  +  l1 - 5 - 

V  4c  /l  3c“  \ 


f  2  4  \r 

b  °"wj  cos  2wt  7  b 2crw  .  0 

—2 — 3M  -  — - rr  sin2' 

v  4c  /  L  6c 


Considering  the  average  torque  per  period  gives  for  the  bracket  contents. 


2t 

-  f 

2t  J 

A  ^  ■ 


d0  =  3/4  ;  0  =  co  t 


so  that 


2  4 

.  -  3tH  b  crco 

AvS  T - 2 - 

16c 


Doing  the  same  for  a  cylinder  of  radius  a  and  subtracting*  for  the  shell 
of  thickness  h  =  (b  -  a)  «  a  gives  the  average  torque  for  length  t 
(Tumbling) 


Avg  T  =  - 
z 


3tH2<7w  t  (b4  -  a4) 


3^H2crwb3f 


♦Since  ^  <p  ^  0  for  the  cylinder,  as  it  was  for  the  sphere,  the  result  now 
is  an  approximation. 
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[4] 

Equation  (113)  agrees  exactly  with  (94)  based  on  Smythe  where  in  (113) 
At  =  1  is  assumed* 


Only  or  *  appears  in  the  average  torque  expression  of  Equation  (107) 
for  a  tumbling  cylinder  shell  and  in  the  torque  expression  of  Equation  (113) 
for  a  spinning  cylinder  shell.  We  now  examine  for  a  possible  case  of 
interest  the  effect  of  using  the  exact  value  of  a 1  as  given  in  the  footnote  on 
page  30  as  opposed  to  the  approximate  value  (for  6  »  b)  also  given  on  the 
bottom  of  page  30  and  used  in  Equations  (107)  and  (113). 


For  a  solid  aluminum  cylinder  with  w  =  2ir/l00  and  of  radius 
b  =  0.625m 


«2'  (approx)  = 


b  cr  co 


0.136 

27 


*2'  (6XaCt)  =  dis¬ 
similarly  for  a  solid  cylinder  of  radius  a  =  0. 615m, 


a  2  (aPProx) 
a  ’  (exact)  = 


0. 132 
2t 


0.127 

27 


The  possibility  of  agreement  appears  plausible  in  the  spinning  case  where 
with  ii|H  a  steady  state  condition  might  result  so  that  the  effect  on  the 

currents  from  V  <p  cancels  motion  effects  because  of  symmetry.  The 
tumbling  case  affords  less  possibility  for  steady  state  conditions. 
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This  amounts  to  a  reduction  of  torque  in  each  case  by  4  percent,  so  that 
using  a  '  (exact)  in  the  shell  cases  gives 

Ld 

revised  Avg  T  for  tumbling  cylinder  shell  =  0. 96  (Avg  T  of 
Equation  (107)) 

revised  T  for  spinning  cylinder  shell  =  0. 96  (T  of  Equation  (113)). 

Finally  with  the  shell  thin  h  «  a,  the  Bessel  function  form  for  a' 
provides  an  approximation  using  the  subtraction  approach  employed  above. 

The  agreement  in  certain  situations  as  we  have  noted,  with  torque  values 
calculated  in  alternate  ways,  e.  g.  using  Smythe  values,  supports  this 
approach,  in  at  least  a  limited  application. 

HEMISPHERE  END  CYLINDER 

Rather  than  open  or  flat  ends  it  may  be  useful  for  radar  applications, 
for  example,  to  have  rounded  ends.  To  provide  some  insight  on  the  effect 
of  such  rounding,  we  choose  to  consider  cylinder  shells  capped  with  hemi¬ 
sphere  shells.  To  relate  the  relative  contribution  of  such  ends  to  the  torque 
we  turn  to  the  dyadic  calculation  in  the  torque  formulation  expressed  by  Tq 
(see  Section  HI). 

a)  Open  End  Cylinder  Shell  Cq 

O 

C  r  r  dV  =  f  Y  x  x.  f.lc. 

J  J  Lj  k  i  k  { 

Vol  k,  i=l 


dV  ;  x  »  x 
x2  =  y 
x„  =  z 


(114) 
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Now 


I  *2  dv  ■  s 


2  ,Tr 

y  dV 


Vol 


Vol 


i/2  b 

C  y  y  (r  cos  0)2  rdfldrdz 
-i/2  a  0 


I 


z  dV  = 

Vol  -i/2  a  0 


iT  f,  4  4*]  . 

=  _  Lb  -  a  J  -  Aq  (115) 

i/2  b  2 7r 

I  j*  S  ^  rd0drdz  =  55  (h2  ”  a2)  =  Bo 


f  x.x 
J  1  j 

Vol 


dV  =  0  ;  all  i  /  j 


(116) 

(117) 


Thus 


I 

,  Vol 


rrdVl  C  =  A  ffx  +  f  i  )  +  B  if 
o  oVxx  y  y  )  ozz 


(118) 
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where 
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Thus 


(i 


rrdV  =  A  I  H  +  H  +  B  H 

J  r  i  V  x  y  1  z 

Vol  /C1  X 


and  as  in  Equation  (120),  we  again  may  write, 


A  (H  +  H  )  +  B,  H  =  A  H.  +  B  Hn 
1  x  yJ  1  z  1  1  1  U 


(125) 


Thus  we  may  use  ^Aq  ,  Bq^  and  ^A^  ,  B^  to  compare  the  relative  effect 
on  torque  of  the  hemisphere  shell  caps. 

We  choose  a  possible  case  of  interest  to  compare  (T)^  and  (T)c  . 

o  1 

For  a  tumbling  cylinder  Cq  or  ,  with  ,  then  over  a  rotation 

period  with  H  fixed,  Avg  T  is  proportional  to  Aq  +  Bq  for  Cq  and 
proportional  to  A^  +  B^  for  since  Tq  is  proportional  to  A^  ( Hj^[ 

+  B^[H|j[^;  k  *  0,  land  [Hjj  *  H  cos  wt,  |H|j|  *  H  sin  ut  . 

CASE  I 

Considering  in  particular  f  =  3m,  b  =  0. 625m  and  h  =  0. 01m  we  get 


A  +  B  =  —  +  B  =  1.25  B  ;  for  C 
oo4  o  o  o 


A  +  B  =  1.25  A  +  2.7B  =  3B  ;  for  C 

11  o  o  o  1 


This  gives  a  factor  increase  of  about  2.4  for  C^  over  C^  . 
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CASE  II 


Again  for  b  =  0. 625m  ,  h  =  0. 01m  and  now  t  =  2. 4m  , 

B 

A  +  B  =  -2  +  B  =  1.33  B  ;  for  C 
o  o  3  o  o  o 


A  +  B  =  1.3  A  +  3.2  B  =  3.6  B  :  for  C 
11  o  o  o  1 


(127) 


This  gives  a  factor  increase  of  about  2. 7  for  C,  over  C  .  Of  course  the 

1  o 

torque  also  is  proportional  directly  to  SL  and  the  overall  length  of  the 
cylinders  is  larger  so  that,  although  the  effect  of  the  ends  being  capped  is 
high  (factor  of  2  to  3  increase  in  torque  for  the  dimensions  used),  a  rough 
comparison  between  cases  1  and  2  gives 


T  (case  2)  2.7  *2  ~ 

T  (case  1)  =  2~i  ^  ~  °'9  (12g. 

Thus,  as  the  length  of  the  cylinder  section  drops,  the  curved  end  caps 
contribute  an  increasing  percentage  to  the  torque.  However,  the  reduction 
in  cylinder  section  length  affects  the  torque  so  as  to  result  in  about  the  same 
torque  for  two  capped  cylinders  of  the  dimensions  considered. 


SECTION  IV 


DYNAMIC  INTERACTIONS  AND  SHELL  DESIGN 


DECAY  IfME  FACTORS 

As  noted  in  the  introduction,  we  are  not  concerned  here  directly  with 
the  time  continuous  motion  of  the  body  interacting  with  the  eddy  current 
torques  which  both  oppose  o>  (decay  torque  component)  and  alter  the 
orientation  of  u?  (turning  torque  component). 

We  concentrate  only  on  the  decay  effects  and  assume  the  u  to  H 
orientation  fixed  in  order  to  assess  in  a  conservative  mamier  the  time 
constants  associated  with  decay.  The  basic  dynamic  equation  relating  decay 
torque  T^  to  acceleration  dw/dt  =  co  is 

TD  =  I  w  (129) 


where  I  is  the  moment  of  inertia  about  the  rotation  axis  in  the  cu  direction. 
We  shall  consider  with  cT 


r]H  (i.e.  H  =  H|) 


D, 


;  k  =  1,  2,  3,  4 


where 

k  =  1:  Spinning  Spherical  Shell 
k  =  2:  Spinning  Cylindrical  Shell  (open-ended) 
k  =  3:  Tumbling  Cylindrical  Shell  (open-ended) 
k  =  4:  Tumbling  Cylindrical  Shell  with  Hemisphere  Shell  Ends 
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For  k  =  1 


From  Equations  (55)  or  (67)  with  3H/3 1  =  0  assumed, 


2t  2  4  - 

— r  H  trb  hw 

3c2 


(130) 


or  more  exactly  with  3H/3 1  =  V  x  E  4  0  ,  from  Equation  (68)  or  (76) 
(i.e.,  at  higher  initial  io  values) 


2  4  — 

2ir  H  trb  hw 


D1  „  T  ,2  2,2  21 

„  2  b  a  h  w 

30  L1  +  - 9 — J 


Equation  (131)  leads  to 


C  OJ 
1 


where 


Then 


dw 

dt  ~  ~~~  o" 

1  +  C.O)^ 
2 


2t  H2  a  b4  h  ,22 ,2  , 

0  =  - - -  ,  c2  =  b  o  h  /9 

3c 


-  +  I  cgj  1  dw  = 
w  2 


-  c  dt 


Taking  w(t  =  0)  =  w  ,  we  get 


£n  («/«o)  +  c2 


/  2  2\ 

/„  -  wo  1  C1 1 


(131) 


(132) 


(133) 


(134) 
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or 


(“/wo) 


(u/“o) 


-  1 


=  e 


-y 


with 


c  .03  c, 

k  =  £  0  k  = _ i 

K1  2  *  2  I 


For  «  1  (i.  e. ,  for  large  wq) 


oj/oj  =  e  y 
o 


So  that  the  exponential  decay  time  constant  becomes 


7  -  i/k  =  1  =  13c 

1  7  2  c,  „  ,,2  ,4 


1  2t  H  a  b  h 


(135) 


(136) 


(137) 


The  same  result  comes  directly  from  Equation  (130)  (the  9H/9 1  =  0  torque 
calculation). 

Thus  it  is  seen  that  t  is  proportional  to  I  and  inversely  proportional 
2  ,  1 

to  H  and  h.  This  will  hold  as  well  for  all  r  .  For  example  with  alumi- 

K 

num  of  spherical  shell  thickness  h  =  ,001m. 

2 

I  =  3.4kg.  m  ,  b  =  0.29m,  H  =  0. 3  oersted*  ,  (see  reference  (7)  )  . 


t  s  70  days. 


*We  take  this  to  be  representative  value  of  earth  magnetic  intensity  at 
about  300  nautical  miles. 
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For  k  =  2 


From  Equations  (94;  or  (113) 

T_  =  -  rr  H2  <j  b3  i  h  w 

D2  - 2 - 

c 

or  more  exactly  from  Equation  (93) 


(138) 


„2  ,  3  .  ,  - 

-  _  rH  ffb  i  hw 

d2' 

c 

Aside  from  new  values  for  c  and  c  ,  Equations  (138)  and  (139)  correspond 
exactly  in  form  to  Equations  (130)  and  (131)  respectively. 

Then  from  Equation  (137) 


1  + 


,2  2  ‘Z  2 
b  cr  h  co 


(139) 


T2 


Ic 


t  H2  a  b3  t  h 


(140) 


For  £  =  3m,  H  =  0.3  oersted,  b  =  0. 625m,  h  =  0. 01m  (1  cm)  and 

2 

aluminum  together  with  an  imposed  value  I  =  75  kg  -  m  (this  value  cannot 
be  attained  with  the  dimensions  given  using  only  aluminum,  see  following 
section). 

s  1.2  days 

We  may  note  in  passing  ihat  for  fixed  H,  and  the  same  material 


(141) 
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For  k  =  3 


From  Equation  (107)  we  obtain  as  an  average  time  constant 


T3  = 


I  16c 


I  4c 


2/4  4\  o  2 

3t  H  (7  ^  b  -  a  )i  3t  H  irb  i  h 

For  the  same  physical  constants  and  dimensions  as  used  in  k  = 


s  1.6  days 


For  k  =  4 


From  Equations  (107)  and  (129)  we  may  write 


Then 


do? 

w 


t  = 
4 


k  =  4 


c^  dt 


k  =  3 


v(°l) 


k  =  4  /  3 


From  Equations  (84),  (126),  and  (12''.'  we  write 


(fj  -  K  (TD) 


k  =  4 


k  =  3 


where,  for  the  dimensions  of  interest, 


K  = 


\  *'k 


=  4 


.  3 


s?  2.5 


(142) 
2  we  get 


(143) 


(144) 

(145) 


(146) 
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For  tumbling  we  consider  I  as  follows  for  k  =  3 


where  y  is  the  mass  density  and  with  volume  element  shown  by 


(147) 


Thus,  using  Equations  (115)  and  (116) 
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« 


In  a  similar  way  but  following  the  type  of  calculation  used  in  (122)  and 
(123),  we  obtain  for  k  =  4 


k  =  4 


so  that 


8;r  (  5  5  \1  7  JE3  (,2 

+  —  l  b  -  a  j  +  - \b  - 

15  '  12  ' 


.y[Al  +  El] 


'A  / Wk . 
SAW*. 


k  =  3  .  1 

-  ,  -  k  E  ’  1 


T4  3  T3 


Therefore  for  the  dimensions  used  in  k  =  3 


=  1. 6  days. 


INERTIA,  WEIGHT  AND  ORBIT  LIFE 

For  purposes  of  shell  design  willi  orbiting  about  the  earth,  we  are 
interested  not  only  in  eddy  current  decay  but  a  host  of  other  effects.  Here 
we  will  concentrate  additionally  only  on  orbit  lifetime  as  determined  by 
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aerodynamic  drag.  *  Also  the  weight  of  the  object  is  of  interest  as  well  as  its 
size.  For  radar  applications  a  conducting  object  is  assumed  and  so  eddy 
current  decays  will  also  be  considered  in  the  design. 


The  model  considered  will  be  as  before  with  the  cylinder  shape  tumbling 
so  that  the  angular  velocity  is  normal  to  H,  but  now  with  the  orbiting  motion 
experiencing  the  effects  of  drag.  The  choice  of  inertia  and  weight  is 
coupled  into  both  the  eddy  current  and  orbit  drag  considerations. 

For  the  purposes  at  hand  we  take 


where 


k  =  3  or  4 


(152) 


“k 

'k 


is  a  constant  depending  on  body  geometry 
is  the  tumbling  moment  of  inertia 


h  is  the  shell  thickness 

Now,  as  noted  previously,  if  I  is  not  arbitrarily  specified 


K3  ’  K4  !4 


(153) 


so  that  only  one  set  of  eddy  current  calculations  need  be  done  for  the  case 
k  =  3  or  k  =  4  cylinders  whenever  I  is  not  imposed. 

Of  course  for  thin  shells,  I  can  be  shown  proportional  to  h.  For 
example,  from  Equation  (149) 


*For  aerodynamic  drag  effects  on  motion  about  the  center  of  mass,  see 
Reference  (19). 
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(154) 


I 

4 


„2  u2 

+  l  b 


+  3£b2 


for  h  «  a.  Thus  t  is  independent  of  h  (with  I  not  imposed). 

A  useful  value  of  r  for  orbit  applications  might  be  about  1.  7  years 
since  this  would  result  in  95  percent  reduction  in  co  in  about  5  years. 

If  we  assume  that  orthogonality  of  H  and  co  is  unlikely,  and  take  H 
2  2 

reduced  to£  0.  8H  (II  to  ^  0.  64  H  )  and  if  we  impose  a  fixed  value  of 
2 

I  =  75  kg-m  ,  then  for  i  =  3m,  b  =  0. 625m,  H  =  0. 3  oersted, 
h  =  0.01m,  as  for  t  ,  we  get  with  k  =  2.4 

O 


say,  for  example 


>  1 

T  =  -  T  3 

4  k  (.  64)  3 


1. 05  days 


=  1.2  days. 

The  thickness  then  required  to  give  =  1.7  years  with  all  other 
factors  remaining  the  same,  becomes 

h.  3  0.0018  cm  3  0.75  mils. 

4 

Now  for  most  operational  radar  applications,  the  radar  skin  depth 


6  = 


2  TOCO 


e 


radary 


t  i  /„  so.  004  cm  =  6  . 

\  1/2  mm, 


so  that 


where  we  desire 


4.5 


6  . 
mm 


h  s  3  6  . 

mm. 
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To  illustrate  a  design,  we  find  using  a  plastic  with  specific  gravity 
=  1.5,  and  h  (aluminum)  =  0.002  cm, 

2 

h  (plastic)  =  1/2  cm:  I  (alum  +  plastic)  =  154  kg-m 

W  (alum  +  plastic)  =  189  lbs. 
t  ~  3.5  years 

2 

h  (plastic)  =  1/4  cm:  I  (alum  +  plastic)  =  76  kg-m 

W  (alum  +  plastic)  =  96  lbs. 

r  s  1.7  years 


Thus  by  striving  to  minimize  torque  by  decreasing  the  conducting 
thickness  while  increasing  the  inertia  with  a  plastic  backing  (without  causing 
a  large  weight  increase)  reasonable  eddy  current  time  constants  for  rotating 
shells  can  be  achieved. 

[  i9 1 

It  may  also  be  mentioned  that  aerodynamic  considerations, 
essentially  independent  of  altitude,  affecting  the  rotating  motion  require 


>  4. 5  in  order  that  the  stable  mode  be  tumbling. 

I  (spin) 

With  the  spin  inertia  calculated  as 

/  4  4  \  o 

I3  (spin)  =  y  ?  \b  -ay  - 

(spin)  »  y  (b4  -  a4)  |  ^  (b5  -  a5)] 

the  above  design  yields 

I  (tumble) 

1  (spin) 


(159) 


(160) 
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5  >  4.5 


Since  we  have  been  varying  weight,  inertia  and  surface  areas,  we  have 
been  changing  the  orbit  lifetime  due  to  aerodynamic  drag. 

Omitting  the  details  of  the  computation  and  considering  eccentric 
orbits  at  300  nautical  miles  perigee  we  find,  using  decay  data,^ 20^  for 
e  =  0  (circular  orbit) 

M 

W  =  200  gives  orbit  life  ss  189  days 

W  =  300^  gives  orbit  life  ss  280  days 

and  for  e  =  0.  05  (apogee  =  694  nautical  miles) 

W  =  200  gives  orbit  life  ss  1658  days 

W  =  300  gives  orbit  life  =£  2430  days 

At  a  circular  orbit  of  450  nautical  miles 

W  =  200  gives  orbit  life  2160  days 

W  =  300  gives  orbit  life  =£  3705  days 

In  these  calculations,  b  =  0.  625  m,  f  =  3m,  an  average  inflight 

2 

projected  area  for  the  hemisphere  capped  cylinder  of  32. 6  ft  and  a 
=  2. 2  were  used. 

It  turns  out  that,  at  the  450-nautical  mile  orbit,  an  orbit  lifetime  about 
twice  as  long  as  the  eddy  current  time  constant  lifetime  is  achieved. 


REFERENCES 


1.  H.  Hertz,  Dissertation  1880,  reprinted  in  Misc.  Papers,  p.  121,  London, 
MacMillan,  1896. 

2.  R.  F.  Gans,  Math  Physik  48,  1  (1903). 

3.  R.  F.  Gans,  Arch  Elektrotech  9,  413  (1921). 

4.  W.  R.  Smythe,  Static  and  Dynamic  Electricity.  2nd  Ed. ,  Chapter  11, 
McGraw  Hill,  1950. 

5.  H.  B.  Rosenstock,  "The  Effect  of  the  Earth's  Magnetic  Field  on  the 
Spin  of  the  Satellite, "  Astronautical  Acta,  3,  pp.  215-221  (1957). 

6.  J.  P.  Vinti,  "Theory  of  the  Spin  of  a  Conducting  Satellite  in  the  Magnetic 
Field  of  the  Earth, "  Ballistic  Missile  Research  Labs,  Aberdeen  Proving 
Ground,  Md. ,  Rept#1020,  July  1957. 

7.  Yu.  V.  Zunov,  "On  the  Problem  of  the  Interaction  Between  a  Satellite 
and  the  Earth's  Magnetic  Field,"  Artificial  Earth. Satellites,  Vol.  3, 
Plenum  Press,  1961. 

8.  R.  P.  Halverson  and  H.  Cohen,  "Torque  on  a  Spinning  Hollow  Sphere 
in  a  Uniform  Magnetic  Field, "  IEEE  Transactions  on  Aerospace  and  . 
Navigational  Electronics,  pp.  IIS-122,  June  1964. 

9.  J.  W.  Hooper,  "The  Damping  of  Metallic  Cylindrical  and  Spherical 
Bodies  Rotating  in  a  Uniform  Magnetic  Field, "  Army  Ballistic  Missile 
Agency,  Huntsville,  Ala.  Tech.  Rept.  DGR-15,  November  1957. 

10.  W.  R.  Smythe,  "Damping  and  Precessing  Torques  on  Spinning  Prolate 
Spheroid  with  Spin  Effect,  "  OOR  Contract  No.  DA-36-034-ORD-1535 
with  Duke  University,  1957. 

11.  R.  H.  Wilson,  Jr. ,  "Geomagnetic  Rotational  Retardation  of  Satellite  1959 
a  1  (Vanguard  HI. 11  Science,  Vol.  131,  Jan  22,  1960. 

12.  R.  J.  Naumann,  "An  Investigation  of  the  Observed  Torques  Acting  on 
Explorer  XI, "  Torques  and  Altitude  Sensing  in  Earth  Satellites, 

Academic  Press,  1964. 


52 


REFERENCES  (Concluded) 


13.  Yu,  E.  Y. ,  "Spin  Decay,  Spin  Precession  Damping,  and  Spin-Axis 
Drift  of  the  Telstar  Satellite, "  Bell  System  Tech.  Jour. ,  pp.  2169- 
2193,  Sept.  1963. 

14.  R.  E.  Robertson,  "Altitude  Control  of  Satellites  and  Space  Vehicles,  " 
Advances  in  Space  Science,  Vol.  2,  Academic  Press,  1960. 

15.  E.  Hecht  and  W.  P.  Manger,  "Magnetic  Altitude  Control  of  the  Tiros 
Satellite. "  Torques  and  Altitude  Sensing  in  Earth  Satellites.  Academic 
Press,  1964. 

16.  L.  D.  Landau  and  E.  M.  Lifshitz,  "Electrodynamics  of  Continuous 
Media, "  Sections  45,  49,  Pergamon  Press,  1960. 

17.  K.  P.  Chopra,  "Interactions  of  Rapidly  Moving  Bodies  in  Terrestial 
Atmosphere. "  Rev,  of  Modern  Physics.  Vol.  33,  No.  2,  pp.  153-189, 
April  1961. 

18.  R.  E.  Fischell,  "Magnetic  Damping  of  the  Angular  Motions  of  Earth 
Satellites, "  ARS  Journal,  pp.  1210-1217,  September  1961. 

19.  K.  R.  Johnson,  "The  Influence  of  Atmospheric  Drag  on  the  Rotational 
Motion  of  a  Satellite,  "  planned  as  a  MITRE  Corp.  MTP. 

20.  J.  B.  Kendrick,  Ed. ,  "Space  Data,  *  TRW  Space  Technology  Labora¬ 
tories,  Redondo  Beach,  California,  1962. 


53 


/ 


Unclassified 

Sri  imtY  Classification 


DOCUMENT  CONTROL  DATA  .RAD 

< ,s»v u/»i)  *  //»*.  m/i«  */ion  •»(  #lile  hotly  oi  abstract  iutd  indexmt J  ii/inofAdun  wii^l  br  entered  *hen  the  overall  reputt  m  t  f/»*wfir-#. 


’  -■>«*  ‘Cltvtlr  fl OrpOfrtK*  Author)  l2«.  REPORT  SCCliniTV  ClA5)'MCAMO'i 

The  MITRE  Corporation  Unclassified 

Bedford,  Massachusetts  »<>■  croup  ~ 


»  Rt'PORt  J  «  H  f. 

EDDY  CURRENT  TORQUES  AND  MOTION  DECAY  ON  ROTATING  SHELLS 


4  r^r  <i.cR«*>r,\,  e  notes  ( Type  of  report  and  tnctusivc  dates) 

N/A 


S  *vThQR.S  ■  t  First  name,  middle  initial,  Jiitnffmf) 

Joseph  F.  A.  Ormsby 


6  REPOPT  OA  *  t 


NOVEMBER  1967 


6m  ~  ON  ’RAf,  t  OP  GRANT  NO 

AF19(628)-5165 

b.  PROJECT  no 


7a.  TOTAL.  NO  OR  PACE5  \7b.  NO  O *  REPS 


9a.  O  RfClN  A  ’'OP’S  REPORT  NuMO£R<5t 


ESD-TR-67-94 


9b.  OTHER  REPORT  NOlSI  CAny  oth«  number',  that  may  be  o'-  ifi.net! 


this  report) 


MTR-224 


JO  O'STPIBUTION  STATEMENT 


This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  unlimited. 


.EMENTARY  notes 


w.  sponsorinc  'Hit**.*  act»v»t>  Directorate  of  Plan- 
ning  and  Technology,  Development  Engineering 
Division,  Electronic  Systems  Division,  L.  G. 
Hanscom  Field,  Bedford,  Massachusetts _ 


’Ttffs  report  describes  a  study  conducted  to  estimate  the  torques  acting  on 
electrically  conducting  shells  rotating  in  a  magnetic  field.  The  basic  electro¬ 
magnetic  expressions  leading  to  an  evaluation  of  torque  are  given  using  various 
approaches  with  comparisons  for  the  spherical  shell.  Application  of  these 
techniques  to  cylindrical  shells  with  flat  and  hemisphere  shell  ends  is  described. 
The  torques  both  slow  and  alter  the  direction  of  the  rotating  motion,  and  the 
resultant  decay  time  constants  can  be  in  the  order  of  days.  Factors  affecting 
this  decay  time  and  influencing  shell  design  are  examined.  ^ 


DD  ,'".“..1473 


Unclassified 

SfiUMlv  CljSSlflfjllliS 


Unclassified 


SeiurtH  I'lasstftruliun 


I  4 

LINK 

A 

UNK  P 

ROLE 

WT 

mm 

W  T 

Eddy  Current 

Electromagnetic  Drag 

Orbit  Decay 

Rotating  Shells 

Satellite  Decay 

m  * 

Unclassified 

Securr.y  Clatsific»u<in 


